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Abstract 

We define higher pentagram maps on polygons in P'^ for any dimension d, which 
extend R. Schwartz's definition of the 2D pentagram map. We prove their integrability 
by presenting Lax representations with a spectral parameter for scale invariant maps. 
The corresponding continuous limit of the pentagram map in dimension d is shown to 
be the (2, d + l)-equation of the KdV hierarchy, generalizing the Boussinesq equation 
in 2D. We also study in detail the 3D case, where we prove integrability for both closed 
and twisted polygons and describe the spectral curve, first integrals, the corresponding 
tori and the motion along them, as well as an invariant symplectic structure. 
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1 Introduction 

The pentagram map was defined by R. Schwartz in [13] on plane convex polygons considered 
modulo projective equivalence. Figure 1 explains the definition: for a polygon P the image 
under the pentagram map is a new polygon T{P) spanned by the "shortest" diagonals of P. 
Iterations of this map on classes of projectively equivalent polygons manifest quasiperiodic 
behaviour, which indicates hidden integrability [H]. 




Figure 1: The image T{P) of a hexagon P under the 2D pentagram map. 

The integrability was proved in [11] for the pentagram map on a larger class of the so 
called twisted polygons in 2D, which are piecewise linear curves with a fixed monodromy 
relating their ends. Closed polygons correspond to the monodromy given by the identity 
transformation. It turned out that there is an invariant Poisson structure for the pentagram 
map and it has sufiiciently many invariant quantities. Moreover, this map turned out to 
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be related to a variety of mathematical domains, including cluster algebras [21 H], frieze 
patterns, and integrable systems of mathematical physics: in particular, its continuous limit 
in 2D is the classical Boussinesq equation [11]. Integrability of the pentagram map for 2D 
closed polygons was established in [13 [12], while a more general framework related to surface 
networks was presented in [3]. 

In this paper we extend the definition of the pentagram map to closed and twisted 
polygons in spaces of any dimension d and prove its various integrability properties. It is 
worth mentioning that the problem of finding integrable higher- dimensional generalizations 
for the pentagram map attracted much attention after the 2D case was treated in [II]0 
The main difficulty in higher dimensions is that diagonals of a polygon are generically skew 
and do not intersect. One can either confine oneself to special polygons (e.g., corrugated 
ones, [3]) to retain the intersection property or one has too many possible choices for using 
hyperplanes as diagonals, where it is difficult to find integrable ones, cf. [9j. 

Below, as an analog of the 2D shortest diagonals for a generic polygon in a projective 
space RP'' we propose to consider a "short-diagonal hyperplane" passing through d vertices 
where every other vertex is taken starting with a given one. Then a new vertex is constructed 
as the intersection of d consecutive diagonal hyperplanes. We repeat this procedure starting 
with the next vertex of the initial polygon. The higher (or d- dimensional) pentagram map 
T takes the initial polygon to the one defined by this set of new vertices. As before, the 
obtained polygon is considered modulo projective equivalence in MP''. 

We also describe general pentagram maps Tp^^ in KP*^ enumerated by two integral param- 
eters p and r by considering p-diagonals (i.e., hyperplanes passing through every pth vertex 
of the polygon) and by taking the intersections of every rth hyperplane like that. There 
is a curious duality between them: the map Tp ^ is equal to modulo a shift in vertex 
indices. However, we are mostly interested in the higher pentagram maps, which correspond 
to T := T2,i in MP'^. 

We start by describing the continuous limit of the higher pentagram map as the evolution 
in the direction of the "envelope" for such a sequence of short-diagonal planes as the number 
of vertices of the polygon tends to infinity. (More precisely, the envelope here is the curve 
whose osculating planes are limits of the short-diagonal planes.) 

Theorem A. (= Theorems 14. 3|, 14.51) The continuous limit of the higher pentagram map 
m m"^ is the {2,d + l)-equation in the KdV hierarchy, which is an infinite- dimensional 
completely integrable system. 

This generalizes the Boussinesq equation as a limit of the pentagram map in MP^ and this 
limit seems to be very robust. Indeed, the same equation appears for an almost arbitrary 

^ There seem to be no natural generalization of the pentagram map to polytopes in higher dimension 
d > 3. Indeed, the initial polytope should be simple for its diagonal hyperplanes to be well defined. In order 
to iterate the pentagram map the dual polytope has to be simple as well. Thus iterations could be defined 
only for d-simplices, which are all projectively equivalent. 
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choice of diagonal planes. It also arises when instead of osculating planes one considers other 
possible definitions of higher pentagram maps (cf. e.g. [9]). 

However, the pentagram map in the above definition with short-diagonal hyperplanes 
exhibits integrability properties not only in the continuous limit, but as a discrete system as 
well. To study them, we define two coordinate systems for twisted polygons in 3D (somewhat 
similar to the ones used in 2D, cf. [Hj), and present explicit formulas for the 3D pentagram 
map using these coordinates (see Theorem 15. 6p . 

Then we describe the pentagram map as a completely integrable discrete dynamical 
system by presenting its Lax form in any dimension and studying in detail the 3D case (see 
Section [6]). For algebraic-geometric integrability we complexify the pentagram map. The 
corresponding 2D case was investigated in [T^ . 

The key ingredient of the algebraic-geometric integrability for a discrete dynamical system 
is a discrete Lax (or zero curvature) equation with a spectral parameter, which in our case 
assumes the following form: 

L,,,+i(A) = P,+i,(A)L,,(A)i^7/(A). 

Here the index t represents the discrete time variable, the index i refers to the vertex of an 
n-gon, and A is a complex spectral parameter. (For the pentagram map in CP'^ the functions 
Lj^t(A) and -Pj,t(A) are matrix-valued of size {d + 1) x [d + 1).) The discrete Lax equation 
arises as a compatibility condition of an over-determined system of equations: 

L,,,(A)^,,i(A) = ^.+i,t(A) 
P,,i(A)*,,i(A) = *,,i+i(A), 

for an auxiliary function \E'j t(A). 

Remark 1.1. Recall that for a smooth dynamical system the Lax form is a differential 
equation of type dfL = [P, L] on a matrix L. Such a form of the equation implies that the 
evolution of L changes it to a similar matrix, thus preserving its eigenvalues. If the matrix 
L depends on a parameter, L = L{X), then the corresponding eigenvalues as functions of 
parameter do not change and in many cases provide sufficiently many first integrals for 
complete integrability of such a system. 

Similarly, an analogue of the Lax form for differential operators of type — L is a zero 
curvature equation dtL — dxP = [P, L] . This is a compatibility condition which provides the 
existence of an auxiliary function ip = ip{t,x) satisfying a system of equations dxip = Lip 
and dfip = Pip . The above Lax form and auxiliary system are discrete versions of the latter. 

In our case, the equivalence of formulas for the pentagram map to the dynamics defined 
by the Lax equation implies complete algebraic-geometric integrability of the system. More 
precisely, the following theorem summarizes several main results on the 3D pentagram map, 
which are obtained by studying its Lax equation. The dynamics is (generically) defined on 
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the space P„ of projectively equivalent twisted ?7,-gons in 3D, which we describe below, and 
has dimension 3n, while closed n-gons form a submanifold of codimension 15 in it. 

Later on we will introduce the notion of spectral data which consists of a Riemann surface, 
called a spectral curve, and a point in the Jacobian (i.e., the complex torus) of this curve, 
as well as a notion of a spectral map between the space Vn and the spectral data. 

Theorem B. (= Theorems 16.151 16.191 17. ip A Zariski open subset of the complexified 
space Vn of twisted n-gons in 3D is a fibration whose fibres are Zariski open subsets of tori. 
These tori are Jacobians of the corresponding spectral curves and are invariant with respect 
to the space pentagram map. Their dimension is 3[n/2\ for odd n and ?){n/2) — 3 for even 
n, where [?t-/2J is the integer part ofn/2. 

The pentagram dynamics on the Jacobians goes along a straight line for odd n and along 
a staircase for even n (i.e., the discrete evolution is either a constant shift on a torus, or its 
square is a constant shift). 

For closed n-gons the tori have dimensions 3[n/2j — 6 for odd n and 3(n/2) — 9 for even 

n. 

Remark 1.2. One also has an explicit description of the the above fibration in terms of 
coordinates on the space of ra-gons. We note that the pentagram dynamics understood as a 
shift on complex tori does not prevent the corresponding orbits on the space Vn from being 
unbounded. The dynamics described above takes place for generic initial data, i.e., for points 
on the Jacobians whose orbits do not intersect certain divisors. Points of generic orbits with 
irrational shifts can return arbitrarily close to such divisors. On the other hand, the inverse 
spectral map is defined outside of these special divisors and may have poles there. Therefore 
the sequences in the space Vn corresponding to such orbits may escape to infinity. 

It is known that the pentagram map in 2D possesses an invariant Poisson structure jTT] . 
which can also be described by using the Krichever- Phong universal formula [TB]. Although 
we do not present an invariant Poisson structure for the pentagram map in 3D, we describe 
its symplectic leaves, as well as the action-angle coordinates. More precisely, we present 
an invariant symplectic structure (i.e., a closed non-degenerate 2- form), and submanifolds 
where it is defined (Theorem 17. 5p . By analogy with the 2D case, it is natural to suggest 
that these submanifolds are symplectic leaves of an invariant Poisson structure, and that the 
inverse of our symplectic structure coincides with the Poisson structure on the leaves. An 
explicit description of this Poisson structure in 3D is still an open problem. 

Note that the algebraic-geometric integrability of the pentagram map implies its Arnold- 
Liouville complete integrability on generic symplectic leaves (in the real case). Namely, the 
existence of a (pre)symplectic structure coming from the Lax form of the pentagram map 
(see Section [7^ . together with the generic set of first integrals, appearing as coefficients of 
the corresponding spectral curve, provides sufficiently many integrals in involution. (Note 
that proving independence of first integrals while remaining within the real setting is often 
more difficult than first proving the algebraic-geometric integrability, which in turn implies 
their independence in the real case.) 



5 



Finally, in Section |8] we present a Lax form for the pentagram maps in arbitrary dimension 
(which implies their complete integrability) assuming their scaling invariance: 

Theorem C. (= Theorem 18. 3p The scale-invariant pentagram map in CP*^ admits a Lax 
representation with a spectral parameter. 

The scaling invariance of the pentagram maps is proved for all d < 6, with some numerical 
evidence for higher values of c? > 6 as well. It would be interesting to establish it in full 
generality. There is a considerable difference between the cases of even and odd dimension 
d, which can be already seen in the analysis of the 2D and 3D cases. 

Acknowledgments. We are grateful to M. Gekhtman and S. Tabachnikov for useful 
discussions. B.K. was partially supported by the Simonyi Fund and an NSERC research 
grant. 

2 Review of the 2D pentagram map 

In this section we recall the main definitions and results in 2D (see [H]), which will be 
important for higher-dimensional generalizations below. We formulate the geometric results 
in the real setting, while the algebraic-geometric ones are presented for the corresponding 
complexification. 

First note that the pentagram map can be extended from closed to twisted polygons. 

Definition 2.1. Given a projective transformation M G P5'L(3,M) of the plane MP^, a 
twisted n-gon in MP^ is a map (j) : Z ^ MP^, such that 0(A; + n) = M o (j){k) for any k. M is 
called the monodromy of (f). Two twisted n-gons are equivalent if there is a transformation 
g G PSL{3, M) such that g o (pi = 02- 

Consider generic n-gons, i.e., those that do not have any three consecutive vertices ly- 
ing on the same line. Denote by Vn the space of generic twisted n-gons considered up to 
PSL{3, M) transformations. The dimension of Vn is 2n. Indeed, a twisted n-gon depends on 
2n, variables representing coordinates of vertices Vk '■= <f){k) for k = 1, ...,n and on 8 param- 
eters of the monodromy matrix M, while the PS'L (3, M)-equivalence reduces the dimension 
by 8. The pentagram map T is generically defined on the space Vn- Namely, for a twisted 
n-gon vertices of its image are the intersections of pairs of consecutive shortest diagonals: 
Tvk ■= (ffc-i, Vk+i) n {vk, Vk+2)- Such intersections are well defined for a generic point in Vn- 

a) Results on integrability (in the twisted and closed cases). There is a Poisson structure 
on Vn invariant with respect to the pentagram map. There are 2[n/2j +2 integrals in 
involution, which provide integrability of the pentagram map on Vn- Its symplectic leaves 
have codimensions 2 or 4 in Vn depending on whether n is odd or even, and the invariant 
tori have dimensions n. — 1 or n. — 2, respectively jllj . 
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Moreover, when restricted to the space C„ of closed polygons (dimCn = 2n — 8 = 2(n— 4)), 
the map is still integrable and has invariant tori of dimension n — 4 for odd n and n — 5 for 
even n. Note that the space C„ of closed polygons is not a Poisson submanifold in the space 
Vn of twisted n-gons, so the corresponding Poisson structure on Vn cannot be restricted to 
C 

There is a Lax representation for the pentagram map. Coefficients of the corresponding 
spectral curve are the first integrals of the dynamics. The pentagram map defines a discrete 
motion on the Jacobian of the spectral curve. This motion is linear or staircase-like depending 
on the parity of n, see |15j . 

b) Coordinates on Vn- The following two systems of coordinates on Vn are particularly 
convenient to work with, see [11]. Assume that n is not divisible by 3. Then there exists 
a unique lift of points Vk = (f){k) G MP^ to the vectors G M.^ satisfying the condition 
det \ Vj, Vj+i, Vj^2\ = 1 for each j. Associate a difference equation to a sequence of vectors 
Vfc G by setting 

Vj+3 = + hjVj+i + Vj 

for all j G Z. The sequences (a-,) and {hj) turn out to be n-periodic, which is a manifestation 
of the fact that the lifts satisfy the relations V^+n = MVj, j G Z, for a certain monodromy 
matrix M G S'L(3, M). The variables a^, < j < n — 1 are coordinates on the space Vn- 
There exists another coordinate system on the space Vn-, which is more geometric. Recall 
that the cross-ratio of 4 points in is given by 

where t is any affine parameter. Now associate to each vertex Vi the following two num- 
bers, which are the cross-ratios of two 4-tuples of points lying on the lines (fj_2,w«-i) and 
{vi+i,Vi+2) respectively: 

Xi = [Vi-2, Vi-i, {{Vi-2,Vi^i) n {Vi,Vi+i)), {{Vi^2,Vi-i) fl {Vi+2, Vi+2))] 

Vi = [((^;^„2,^^^-l) n {vi+i,Vi+2)), {{vi-i,Vi) fl {vi+i,Vi+2)), Vi+i, Vi+2] 
In these coordinates the pentagram map has the form 

1 - Xi-iVi-i 1 - Xi+2yi+2 

T Xi = Xi- T yi = Vi+i— . 

1 - Xi+iUi+i 1 - XiUi 

One can see that the pentagram map commutes with the scaling transformation pTl [T3] : 

Rs ■ {Xl,yi, Xn, Vn) {SXI, S~^yi, SXn, S'^^Un) - 

In these coordinates the invariant Poisson structure has a particularly simple form, see [TT] . 
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c) Continuous limit: the Boussinesq equation. The n ^ oo continuous limit of a twisted 
n-gon with a fixed monodromy M G PSL{3,M.) can be viewed as a smooth parameterized 
curve 7 : M — )■ MP^ satisfying 7(0; + 2n) = M'~f{x) for all a; G M. The genericity assumption 
that every three consecutive points of an n-gon are in general position corresponds to the 
assumption that 7 is a non-degenerate curve in MP^, i.e., the vectors 'j'{x) and 'j"{x) are 
linearly independent for all x G M. 



G{x) 




Figure 2: Constructing the envelope L^{x) in 2D. 



The space of such projectively equivalent curves is in one-to-one correspondence with lin- 
ear differential operators of the third order: L = + Ui{x)d + Uq{x), where the coefficients 
Uq and Ml are periodic in x. Namely, a curve 7(0;) in MP^ can be lifted to a quasi-periodic 
curve G = {G{x)} in satisfying det \G{x),G'{x),G"{x)\ = 1 for all x G M. The compo- 
nents of the vector function G{x) = {Gi{x),G2{x),G3{x)) are homogenous coordinates of 
7(x) in RP^: j{x) = (Gi : G2 : G3){x) G MP^ The vector function G{x) can be identified 
with a solution of the unique linear differential operator L, i.e., the components of G{x) are 
identified with three linearly independent solutions of the differential equation Ly = 0. 

A continuous analog of the pentagram map is obtained by the following construction. 
Given a non-degenerate curve 7(x), we draw the chord (7(0; — e),'j{x + e)) at each point 
X. Consider the envelope ie{x) of these chords. (Figure 2 shows their lifts: chords {G{x — 
e),G{x + e)) and their envelope L^{x).) Let ui^e and Mo,e be the periodic coefficients of the 
corresponding differential operator. Their expansions in e have the form Ui^e = Ui + e^Wi + 
C(e^) and allow one to define the evolution dui/dt := Wi, i = 0,1. After getting rid of Uq 
this becomes the classical Boussinesq equation on the periodic function u = Ui, which is the 
(2, 3)-fiow in the KdV hierarchy of integrable equations on the circle: Utt+'2{u'^)xx+Uxxxx = 0. 

Below we generalize these results to higher dimensions. 
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3 Geometric definition of higher pentagram maps 



3.1 Pentagram map in 3D 

First we extend the notion of a closed polygon to a twisted one, similar to the 2D case. We 
present the 3D case first, before giving the definition of the pentagram map in arbitrary 
dimension, since it is used in many formulas below. 

Definition 3.1. A twisted n-gon in MP^ with a monodromy M G >S'L(4, M) is a map : Z — )■ 
MP'^, such that + n) = M o (^{k) for each /c G Z. (Here we consider the natural action of 
5'L(4,M) on the corresponding projective space MP^.) Two twisted n-gons are (projectively) 
equivalent if there is a transformation g G >S'L(4, M), such that g o (pi = 02- 

Note that equivalent n-gons must have similar monodromies. Closed n-gons (space poly- 
gons) correspond to the monodromies M = Id and —Id. Let us assume that vertices of an 
n-gon are in general position, i.e., no four consecutive vertices belong to one and the same 
plane in MP^ Also, assume that n is odd. Then one can show (see Section 15.31 below and 
cf. Proposition 4.1 in [llj) that there exists a unique lift of the vertices Vk ■= (j){k) G MP^ 
to the vectors Vk G satisfying for all j G Z the identities det \ Vj, V^+i, Vj+2, V^+sl = 1 and 
Vj^n = MVj, where M G 5'L(4, M)|§ These vectors satisfy difference equations 

Vj+i = ajVj+3 + bjVj+2 + CjVj+i - Vj, j G Z, 

with n-periodic coefficients {aj, bj, cj) and we employ this equation to introduce the (a, b, c)- 
coordinates on the space of twisted ra-gons. 

Define the pentagram map on the classes of equivalent n-gons in such a way. 

Definition 3.2. Given an n-gon in RP^, for each G Z consider the two-dimensional 
"short-diagonal plane" Pk := {vk-2,Vk,Vk+2) passing through 3 vertices ffc_2, ffc, ffc+2- Take 
the intersection point of the three consecutive planes Pk-i, Pk, Pk+i and call it the image 
of the vertex Vk under the space pentagram map T, see Figure 3. (We assume the general 
position, so that every three consecutive planes Pk for the given n-gon intersect at a point.) 

By lifting a two-dimensional plane Pk from MP^ to the three-dimensional plane through 
the origin in (and slightly abusing notation) we have Pk = *{Vk-2 A A Vk+2) in terms 
of the natural duality * between R'^ and R'^*. The lift of Tvk to is proportional to 

*[Pk-i A Pk A Pk+i]. 

Below we describe the properties of this space pentagram map in detail. 

^This explains our choice of the group 5i(4, R) rather than the seemingly more natural group PSL{A, R): 
since 5i(4,R) is a two-fold cover of P5'L(4,M), a twisted n-gon would have two different lifts from MP'^ to 
R"* corresponding to two different lifts of the monodromy M from the latter group. 
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Vk 

Figure 3: In 3D the image Tvk of the vertex Vk is the intersection of three "short-diagonal" planes 
Pk-i,Pk, and Pk+i- 

3.2 Pentagram map in any dimension 

Before defining the pentagram map in MP"^, recall that SL{d + 1, M) is a two- fold cover of 
PSL{d + 1, M) for odd d and coincides with the latter for even d. 

Definition 3.3. A twisted n-gon in RP'^ with a monodromy M G SL{d + 1,M) is a map 
: Z — 7- MP'^, such that (j){k + n) = M o (f){k) for each A; G Z, and where M acts naturally on 
MP''. 

We define the SL{d + 1, M)-equivalence of ra-gons as above, and assume the vertices 
Vk '■= 4>{k) to be in general position, i.e., in particular, no d + 1 consecutive vertices of an 
n-gon belong to one and the same {d — l)-dimensional plane in MP*^. 

Remark 3.4. One can show that there exists a unique lift of the vertices Vk = 4>{k) G MP'' 
to the vectors Vk G M''"'"^ satisfying det |Vj-, Ij+i, V^+dl = 1 and = MVj, j G Z, where 
M G SL{d + 1, M), if and only if the condition gcd{n, d + 1) = 1 holds. The corresponding 
difference equations have the form 

Vj+d+i = o-j^'^j+d + aj^d-iVj+d-i + ••• + a.j,iVj+i + {-lYVj, j G Z, (1) 

with n-periodic coefficients in the index j. This allows one to introduce coordinates {aj^k, < 
j < n — 1, 1 < k < d} on the space of twisted n-gons in MP'^. 

For a generic twisted n-gon in MP'^ one can define the "short-diagonal" (c?— l)-dimensional 
plane Pk passing through d vertices of the n-gon by taking every other vertex starting at 
the point Vk, i.e., through the vertices Vk,Vk+2, ■■■,Vk+2d-2- For calculations, however, it is 
convenient to have the set of vertices "centered" at Vk, and then the definition becomes 
slightly different in the odd and even dimensional cases. 
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Namely, for odd dimension d = 2x + 1 we consider the short- diagonal hyperplane Pk 
through the d vertices 

(thus including the vertex Vk itself), while for even dimension d = 2x we take Pk passing 
through the d vertices 

Pk '■= {Vk~2K+l,Vk~2K+3, ■■■,Vk-l,Vk+l, ■■■,Vk+2K-l) 

(thus excluding the vertex Vk). 

Definition 3.5. The higher pentagram map T takes a vertex Vk of a generic twisted n-gon 
in MP*^ to the intersection point of the d consecutive short-diagonal planes Pi around Vk- 
Namely, for odd d = 2x + 1 one takes the intersection of the planes 

Tvk := Pk-^ n Pk-^+i n ... n Pfe n ... n Pk+^ , 

while for even d = 2x one takes the intersection of the planes 

Tvk := Pk-^+i n Pk-^+2 n ... n Pfc n ... n Pk+^ . 

The corresponding map T is well defined on the equivalence classes of n-gons in MP*^. 

As usual, we invoke the generality assumption to guarantee that every d consecutive 
hyperplanes Pj intersect at one point in MP'*. It turns out that the pentagram map defined 
this way has a special scaling invariance, which allows one to prove its integrability: 

Theorem C. (= Theorem 18. 3|) The scale-invariant higher pentagram map is a discrete 
completely integrable system on equivalence classes ofn-gons in MP"'. It has an explicit Lax 
representation with a spectral parameter. 

Remark 3.6. The pentagram map defined this way in ID is the identity map. In the 2D 
case this definition was given in [13] and its integrability for twisted polygons was proved in 
[TT] , while for closed ones it was proved in [151 [12] • 

Remark 3.7. One can also give an "asymmetric definition" for planes Pk, where more 
general sequences of d vertices Vkj are used, and then Tvk is defined as the intersection of 
d consecutive planes Pk- It turns out, however, that exactly this "uniform" definition of 
diagonal planes Pk, where Pk passes through every other vertex, leads to integrability of the 
pentagram map. 

One of possible definitions of the pentagram map discussed in [H] coincides with ours in 
3D. In that definition one takes the intersection of a (possibly asymmetric, but containing 
the vertex Vk) plane Pk and the segment [vk-i,Vk+i]: for our centered choice of Pk this 
segment belongs to both planes Pk-i and Pk+i- The definitions become different in higher 
dimensions. 
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3.3 General pentagram maps and duality 



We define more general pentagram maps Tp^r depending on two integral parameters in ar- 
bitrary dimension d. These parameters specify the diagonal planes and which of them to 
intersect. 

Definition 3.8. For a generic twisted n-gon in MP*^ one can define a p-diagonal hyperplane 
Pk as the one passing through d vertices of the n-gon by taking every pth vertex starting at 
the point Vk, i.e., 

Pk '■= {Vk, Vk+p, Vk+{d-l)p) ■ 

The image of the vertex Vk under the general pentagram map Tp^^ is defined by intersecting 
every rth out of the p-diagonal hyperplanes starting with Pk. 

Tp^rVk '■= Pk n Pk+r n ... n Pk+{d-l)r ■ 

The corresponding map Tp^r is considered on the space of equivalence classes of n-gons in 
MP''. 

For the higher pentagram map T discussed in Section one has p = 2, r = 1, and the 
indices in the definition of Pk are "centered" at Vk- In other words, T = T2^i o Sh., where 
Sh stands for some shift in the vertex index. Below we denote by Sh any shift in the index 
without specifying the shift parameter. Note that Tp^p = Sh. 

Theorem 3.9. There is a duality between the general pentagram maps Tp^r and T^^p: 

Tp^r = T~p oSh. 

For example, the map Ti^2 in 2D is defined by extending the sides of a polygon and 
intersecting them with the "second neighbouring" sides. This corresponds exactly to passing 
from T{P) back to P in Figure 1, i.e., it is the inverse of T modulo the numeration of vertices. 

Proof. To prove this theorem we introduce the following "duality maps," cf. |TT]. 

Definition 3.10. Given a generic sequence of points 0(j) G MP"', j E Z, and a nonzero 
integer p we define the sequence ap{(j){i)) G (MP'^)* as the plane 

otpW)) := <PU + P), 0(J + {d- l)p)) . 

The following proposition is straightforward for our definition of the general pentagram 
map. 

Proposition 3.11. For every nonzero p the maps ap are involutions modulo index shifts 
(i.e., ttp = Sh), they commute with the shifts (i.e., ap o Sh = Sh o ap), and the general 
pentagram map is the following composition: Tp^r = a^ o apO Sh. 
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To complete the proof of the theorem we note that 

T~p = (ttp o o Sh)~^ = Sh~^ oa~^ oa~^ = ShoarO ShoapO Sh = arOapO Sh = Tp^r ° Sh, 

since = OpO Sh from the Proposition above, while Sh"^ = Sh and Sh o Sh = Sh. □ 

Below we construct a Lax form for the higher pentagram maps, i.e., for the maps T2^i 
(and hence for Ti 2 as well) for any d. Integrability of the pentagram map on a special class 
of the so-called corrugated twisted polygons in MP'^ was proved in [3j, which should imply 
the integrability of the pentagram map Tp^i in 2D. Then the above duality would also give 
integrability of Ti p in MP^, defined as the intersection of a pair of polygon edges whose 
numbers differ by p. (One should also mention that for p and r mutually prime with n one 
can get rid of one of the parameters by appropriately renumbering vertices, at least for the 
closed ra-gon case, cf. [12] in 2D. This reduces the study to the map Tp^i for some values of 
n.) Complete integrability of general pentagram maps for other pairs {p, r) in WP'^ is a wide 
open problem. 

Problem 3.12. Which of the general pentagram maps Tp^r I^P'^ ore completely integrable 
systems? 

4 Continuous limit of the higher pentagram maps 
4.1 Definition of the continuous limit 

In this section we consider the continuous limit of polygons and the pentagram map on them. 
In the limit n — )■ cxo a twisted ra-gon becomes a smooth quasi-periodic curve 7(x) in MP*^. 
Its lift G{x) to M^+^ is defined by the conditions: i) the components of the vector function 
G{x) := {Gi{x), ...,Gd+i{x)) provide homogeneous coordinates for 7(x) = (Gi : ... : Gd+i){x) 
in RP"^, ii) det \G, G', G^'^^Kx) = 1 for all x G M, and Hi) G{x + 2tt) = MG{x) for a given 
M G SL{d + 1, M). Then G satisfies a linear differential equation of order d + 1: 

G^'^+^^ + Ud-i{x)G^'^-^^ + ... + ui{x)G' + uo{x)G = (2) 

with periodic coefficients Ui{x). Here and below ' stands for d/dx. 

Let us consider the case of odd d = 2x+ 1. Fix a small e > 0. A continuous analog of the 
hyperplane Pk is the hyperplane Pe{x) passing through d points 7(0; — xe), 7(x), 7(x + 
xe) of the curve 7^1 Note that as e — )■ the hyperplanes Pe{x) tend to the osculating 
hyperplane of the curve 7 spanned by the vectors 7'(a;), 7"(x), 7'^'^~^^(a;) at the point 7(0;). 

Let ie{x) be the envelope curve for the family of hyperplanes Pe{x) for a fixed e. The enve- 
lope condition means that for each x the point ie{x) and the derivative vectors 

For a complete analogy with the discrete case, one could take the points j{x — 2>fe), 
7(0; — 2{k — l)e), 7(0;), 7(x + 2k€). However, one can absorb the factor 2 by rescaling e — >■ 2e. 
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belong to the plane Pe{x). This means that the lift of ie{x) to L^{x) in M*^"*"^ satisfies the 
system of c? = 2x + 1 equations (see Figure 4): 

det\G{x - >ce),G{x - {x - l)e), ...,G{x), ...,G{x + xe),Li^\x)\ =0, j = 0, ci - 1. (3) 



G{x) 



G{x-e) 




Figure 4: The envelope L^{x) in 3D. The point L^{x) and the vectors L'i,{x) and L'l{x) belong to 
the plane {G{x - e), G(x + e)). 

Similarly, for even d = 2x the lift L^{x) satisfies the system of d equations: 
det 1^(3: - (2x - l)e), G{x - (2x - 3)e), G{x - e), G{x + e), ... 

...,G(x + (2x-l)e),L(^)(x)| =0, j = 0,...,d-l. (4) 

The evolution of the curve 7 in the direction of the envelope as e changes, defines a 
continuous limit of the pentagram map T. Namely, below we show that the expansion of 
L^{x) has the form 

L,{x) = G{x) + e^B{x) + C(e^). 
The family of functions L^{x) satisfies a family of differential equations: 

L^'^+i) + + ... + u,,,{x)L[ + Mo,.(x)L, = 0, 

where Ujfl{x) = Uj{x). 

Expanding the coefficients Uj^^ix) as Uj^e{x) = Uj{x) + e^Wj{x) + C(e^), we define the 
continuous limit of the pentagram map T as the system of evolution differential equations 
duj{x)/dt = Wj{x) for j = 0, d — 1, i.e., plays the role of time. 

Theorem A. The continuous limit duj{x)/dt = Wj{x), j = 0,...,d — 1 for x E of the 
pentagram map is the {2,d+ l)-KdV equation of the Adler-G elf and- Dickey hierarchy on the 
circle. 

This theorem is proved as a combination of Theorems 14.31 and 14.51 below. 
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Remark 4.1. Recall the definition of the KdV hierarchy (after Adler-Gelfand-Dickey, fL\). 
Let L be a linear differential operator of order d + 1: 

L = S'^+i + ua-iix)d''-' + Ud^2ix)d''-^ + ... + u,ix)d + uoix) (5) 

with periodic coefficients Uj{x), where (9^ stands for d^/dx'^. One can define its fractional 

power 1/^/'^+^ as a pseudo-differential operator for any positive integer m and take its purely 

2 

differential part ■= . In particular, for m = 2 one has Q2 = + — 

Then the (m, d+ l)-KdV equation is the following evolution equation on (the coefficients of) 
L: ^ 

-L=[Q^,L]. 

Remark 4.2. a) For d = 1 the discrete pentagram map is the identity map, hence the 
continuous limit is trivial, which is consistent with vanishing of the (2,2)-KdV equation. 

b) For d = 2 the (2,3)-KdV equation is the classical Boussinesq equation, found in [TT] . 

c) Apparently, the (2, (i+ l)-KdV equation is a very robust continuous limit. One obtains 
it not only for the pentagram map defined by taking every other vertex, but also for a non- 
symmetric choice of vertices for the plane Pk, see Remark 14.41 Also, the limit remains the 
same if instead of taking the envelopes one starts with a map defined by taking intersections 
of various planes |9j. 



4.2 Envelopes and the KdV hierarchy 

Theorem 4.3. For any dimension d, the envelope L^{x) has the expansion 

L,{x) = G{x) + Cd {g'\x) + -^Ud-iG{x)^ + O(e^) 

for a certain constant Cd, as e — )■ 0. 

The e^-term of this expansion can be rewritten as Cd + ^ ^ ^ Ud-i{x)^ G{x). Conse- 
quently, it defines the following evolution of the curve G{x): 



Proof. Since approaches G as e — > 0, one has the expansion = G + eA + e^B + 
e^C + 0{e^). First we note that the expansion of in e has only even powers of e, since the 
equations ([3]) and (jl]) defining have the symmetry e — — e. Therefore, we have A = G = 
and L, = G + e^B + 0{e^). 
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Notice that G{x) with its first d derivatives form a basis in M'^+^ for each x. We express 
the vector coefficient B in this basis: B = bQ{x)G + hi{x)G' + ... + bd{x)G^'^\ Recall that, 
e.g., for odd d = 2x + 1 the lift L^{x) satisfies the system of d equations: 

det\G{x- xe),G{x-{>c-l)e),...,G{x),...,G{x + >ce),Li^\x) | = 0, j = 0,...,d-l. 

2 

Fix X and expand all terms in e: e.g., G{x + e) = G{x) + eG'{x) + yC" + etc. In each 
equation consider the coefficients at the lowest power of e, being 2 + d{d — l)/2 here. 

The equation with j = gives Gd det \G, G', G^'^~^\ B\ = for some nonzero Gd, which 
implies that there is no G'^'^^-term in the expansion of B. Similarly, for j = 1 we obtain 
Crfdet \G, G', G^'''~^\ B'\ = 0, which means that there is no G'^'^^-term in the expansion of 
B', or, equivalently, there is no -term in the expansion of B. Using this argument for 

j = 0, 1, ...,d — 3, and c? — 1, we deduce that B contains no terms with G^'^\ G'" and G'. 

The equation with j = d — 2 results in a different term in the expansion and gives 

det \G, G', G^'^-'\ 5(^^-2) I _ Gd det \G, G\ G^'^-'\ G^^^ | , 

which implies that B = GdG" + b{x)G for some function b{x). 

Finally, the normalization det | L^, L^, Li"^^ \ = 1 allows one to find b{x) by plugging in 
it = G + e^{GdG" + b{x)G) + C(e^). For the e^-terms one obtains 

Gd det I G, G', G^''~^\ G^'^+i), G^''^ \ + Gd det | G, G', G^'^-^), G^''~'\ G('^+2) | 

+{d + l)b{x) det I G, G', G('^-2), G(^-i), G^'*) | = . 

By using the linear differential equation G^'^'^^^ + Ud-i{x)G^'^~^^ + ... + Uq{x)G = to express 
Q{d+i) Q^y^^ Q{d+2) yjg^ lower derivatives we see that the first and the second determinants are 
equal to — m^-i, while the last one is equal to 1. Thus one has {d + l)b{x) — 2Gd Ud-i{x) = 0, 
which gives b{x) = Gd-^Ud-iix). 

Hence = G + G^ (G" + ^Ud^iG) + 0{e^), as required. □ 



Remark 4.4. One can see that the only condition on vertices defining the hyperplane Pk 
required for the proof above is that they are distinct. A different choice of vertices for the 
hyperplane Pk changes the constant Gd, but does not affect the evolution equation for G. 
The above theorem for an envelope L^{x) is similar to an analogous expansion in [9j for a 
curve defined via certain plane intersections. 

Theorem 4.5. In any dimension d the continuous limit of the pentagram map defined by 
the evolution 

of the curve G coincides with the {2,d+ l)-KdV equation. Consequently, it is an infinite- 
dimensional completely integrable system. 
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Proof. Recall that the (2, d + l)-KdV equation is defined as the evolution 



L=[Q2,L] := Q2L-LQ2, 

where the linear differential operator L of order c? + 1 is given by formula ([5]) and Q2 = 
2 

+ Mrf-i- Here L stands for dL/dt. 

d + l ' 
By assumption, the evolution of the curve G is described by the differential equation 

G = Q2G. We would like to find the evolution of the operator L tracing it. For any t, the 

curve G and the operator L are related by the differential equation LG = of order d + l. 

Consequently, 

LG + LG = 0. 

Note that if the operator L satisfies the (2, d + l)-KdV equation and G satisfies G = Q2G, 
we have the identity: 

LG + LG= {Q2L - LQ2)G + LQ2G = Q2LG = . 

In virtue of the uniqueness of the linear differential operator of the form (|5]) for a given 
fundamental set of solutions, we obtain that indeed the evolution of L is described by the 
(2, d + l)-KdV equation. □ 



Remark 4.6. The proof above is reminiscent of the one used in (TU] to study symplectic 
leaves of the Gelfand-Dickey brackets. Note that the absence of the term linear in e is related 
to the symmetric choice of vertices for the hyperplane Pk- For a non-symmetric choice the 
evolution would be defined by the linear term in e and given by the equation G = G'. This 
is the initial, {l,d + l)-equation of the corresponding KdV hierarchy, manifesting the fact 
that the space x-variable can be regarded as the "first time" variable. A natural question 
arises whether the whole KdV-hierarchy is hidden as an appropriate limit of the pentagram 
map. An evidence to this is given by noticing that the terms with higher powers in e lead 
to equations similar to the higher equations in the KdV hierarchy, see Appendix 19. 2[ 

Remark 4.7. One can see that the continuous limit of the general pentagram maps T^ ,. 
for various p ^ r in RF'^ defined via envelopes for a centered choice of vertices is the same 
(2, d + l)-KdV fiow, i.e., the limit depends only on the dimension. 

Indeed, an analog of the p-diagonal is the plane Pe(x) passing through the points 
G{x),G{x + ep), G{x + e{d — l)p). Rescaling e, we can assume the points to be G{x), G{x + 
e), G{x + e{d— 1)), which leads to the planes Pe{x) defined in Section l^?T] after a shift in x. 
Then the definition of L^{x) via the envelope of such planes will give the same (2, d+l)-KdV 
equation. 

It would be interesting to define the limit of the intersections of every rth plane via some 
higher-order terms of the envelope, as mentioned in the above remark, so that it could lead 
to other (m, d + l)-equations in the KdV hierarchy. 
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5 Explicit formulas for the 3D pentagram map 
5.1 Two involutions 

Now we return to the 3D case. In this section we assume that n is odd and consider n-gons in 
MP^ Recall that in this case an n-gon with a given monodromy M G 5'L(4, M) lifts uniquely 
to and is described by difference equations 

Vj+i = ajVj+3 + bjVj+2 + CjVj+i - Vj, j G Z, (6) 

with n-periodic coefficients (a^, bj, Cj). In other words, for odd n the variables {aj, bj, cj), < 
j < n — 1, provide coordinates on the space Vn of twisted n-gons in RP^ considered up to 
projective equivalence (see Proposition 15. 9p . 

In order to find explicit formulas for the pentagram map, we present it as a composition 
of two involutions a and /3, cf. Section 13.31 and then find the formulas for each of them 
separately. The same approach was used in [TT] in 2D, although the formulas in 3D are more 
complicated. 

Definition 5.1. Given a sequence of points G MP^, j G Z, define two sequences 

e (P3)* and e (MP^)*, where 

«) is the plane (0(j - l),0(j),0(j + 1)); 

b) is the plane (0((j - 2), <^(i), 0(j + 2)). 

Proposition 5.2. The maps a and /3 are involutions, i.e., = = Id, while the pentagram 
map is their composition: T = a o /3. 

Note that the indices which define the planes are symmetric with respect to j. As a 
result, we do not have an extra shift of indices, cf. Proposition 13.111 (unlike the 2D case and 
the general map ap). 

Lemma 5.3. The involution a : Vi ^ Wi = *{yi /\ Vi^i A Vj+i) maps equation ^ to the 
following difference equation: 

Wi+^ = Ci+iWi+3 + biWi+2 + ai-iWi+i - Wi. 

Lemma 5.4. The involution (3 : Vi ^ Wi = \i * {Vi /\ Vj_2 A rnaps equation ^ to the 

difference equation 

where the coefficients Ai, Bi, Ci are defined as follows: 

Ai = Cj_i(ajaj4.2 + aj+2&t+iCj + K+iK+2^'i+4:, 
Bi = {{ai-2 + 6j-iCi_2)(ci+2 + ai+2bi+i) - ai+2Q_2)x 
X (ai-iCj+i + di+ibiCi-i + Cj_iCj+i)AjAi+iAj+3Aj+4, 
Ci = aj+i(ajaj_(.2 + ai+2&i+iCj + CjCj+2)^A.j4_2Aj+3A^_,_4. 
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The sequence Aj, i ^ is n-periodic and is uniquely determined by the condition 



(ai_2ai + aibi^iCi-2 + Q_2Q)(ai_iai+i + ai+ibid-i + Ci-iCi+ij 

The proofs of these lemmas are straightforward computations, which we omit. Combined 
together, these lemmas provide formulas for the pentagram map. They have, however, one 
drawback: one needs to solve a system of equations in Aj, i G Z, which results in the 
non-local character of the formulas in (a, b, c)-coordinates and their extreme complexity. 



5.2 Cross-ratio type coordinates 

Similarly to the 2D case, there exist alternative coordinates on the space Vn- They are 
defined for any n, and the formulas for the pentagram map become local, i.e., involving the 
vertex itself and several neighboring ones. 

Definition 5.5. For odd n the variables 

= — , Vj = , Zj = , < j < - 1, 

Ojj^\Cj (Xjj^\Oj 

provide coordinates on the space where the n-periodic variables {aj,bj,Cj), j G Z, are 
defined by the difference equation ([6]). 

It turns out that the variables {xj, yj, Zj), < j < n — 1 are well defined and independent 
for any n, even or odd. Below we provide two (equivalent) ways to define them for even 
n: a pure geometric (local) definition of these variables (see Proposition 15.71) and the above 
definition extended to quasi-periodic sequences {aj,bj,Cj) (see Section [5731) . 

Theorem 5.6. In the coordinates Xi,yi,Zi the pentagram map for any (either odd or even) 
n is given by the formulas: 



rj.*i N _ 1 + Vi-l + + yi-lZi+2 - Vi+lZ, 



T\y^ 
T*izi 



1 + Vi-i + Zi 

Xi_iyi_iZi (1 + yj+i + z»+2)(l + yi-2 + Zi-i) 

XiZi-i {l + yi + Zi+i){l + yi_i + Zi+2 + yi-iZi+2 - yi+iZi^ 

Xi+iZj (1 + yi+i + 2:^+2) (1 + yi-2 + Zi_i) 

Xi (1 + yi_i + Zi){l + yi_2 + Zi+i - yiZi^i + yi-2Zi+i) ' 



Before proving this theorem we describe the {xi,yi,Zi) coordinates in greater detail. It 
turns out that they may be defined completely independently of {ai,bi,Ci) in the following 
geometric way. 

Recall that the x, y coordinates for the 2D pentagram map are defined as cross-ratios for 
quadruples of points on the line {Vi, Vi+i), where two points are these vertices themselves, and 
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two others are intersections of this hne with extensions of the neighbouring edges. Similarly, 
the next proposition describes the new coordinates via cross-ratios of quadruples of points, 
2 of which are the vertices of an edge, and 2 others are the intersection of the edge extension 
with two planes. For instance, the variable yi is the cross-ratio of 4 points on the line 
{Vi, l^+i), two of which are Vi and l^+i, while two more are constructed as intersections of 
this line with the planes via the triple (Vi+2, V^+s, Vi+4) and with the plane via the triple 
(1^+2, Vi+4, Vi+5). More precisely, the following proposition holds. 

Proposition 5.7. The coordinates Xi,yi, Zi are given by the cross-ratios: 

= — V'j+s, $012 5 ^123]^ 
y, = -[V,V+^,<^ll„<^ll,], 

z, = -[V+,,V+,,<!>t>n,^?23l 

where the point $:^f^m2,m3 /^'^ given i is the intersection of the line (Vi+j^, V^+jJ '^'^^^ 
plane (V^+mi, ^^+^2? Vi^ms)- 

By the very definition these coordinates are projectively invariant. 



Proof of proposition. If * is the Hodge star operator with respect to the Euclidean metric 
in M^(or C^), then 

*mf'm2,m3 = * A Vi+,,) A A V+r^, A ^^+^3)) . 

It suffices to prove the proposition in the case of an odd n, because the formulas are local, 
and we can always add a vertex to change the parity of n. Therefore, we may assume that 
{ttj, bj, Cj) are global coordinates and use them for the proof. 
A simple computation shows that 

*oi2 = -{bi+i + aiai+i)Vi+4 + aiVi+5, 

*234 = V^- C,Vi+i, 

*245 = -bi+iV + (ai + Cibi+i)Vi+i, 
*oi3 = + biai+i)Vi+4. - biVi+5, 

$123 = —(^i+lVi+i + 

Recall (see Lemma 4.5 in 111]) that if 4 vectors a,b,c, d e M.^ (or C^) lie in the same z- 



dimensional plane and are such that 

c = Xia + A2&, d = fiia + fi2b, 
then the cross-ratio of the lines spanned by these vectors in the plane is given by 

[a,6,c,rfj = . 

-^2/^1 



20 



Comparing the cross-ratios with the original definition of the variables Xi,yi^Zi concludes 
the proof. □ 

Now we are in a position to prove the explicit local formulas in Theorem 15.61 The proof 
is similar, but more involved than that of Proposition 4.11 in [TT] . 

Proof of theorem. Due to the local character of the formulas for the pentagram map 
T in (xj, 2;j)-coordinates, we may always add an extra vertex to make the number n of 
vertices odd, and then use coordinates {aj,bj,Cj) and Lemmas 15.31 and for the proof. 
The pentagram map is a composition T = a o /3 : — )■ f/j. Namely, 

U^ = fii* [*{Vi^3 A A V^+i) A *{V,-2 AV^A Vi+2) A *{Vi.i A V^+i A Vi+3)] , 

where the constants /Xj are chosen so that det \ Uj, Uj+i, f^j+sl = 1 for all j. 
At the level of the coordinates (a^, bj, Cj), we obtain: 



T [Xi) = - — - — , T [yi) = — — , T [z, 



Ci+lCi+2 -Bi+iv4j_i Ci^2Bi 



where Ai,Bi,Ci are defined in Lemma [5.4[ Eliminating the variables Aj with different i by 
using the formula for the product AjAj+iAj+2Aj+3 concludes the proof. □ 



5.3 Coordinates on twisted polygons: odd vs. even n 

In this section we compare how one introduces the coordinates on the space Vn of twisted 
n-gons for odd or even n, and how this changes the statements above. 

Definition 5.8. Call a sequence {aj,bj,Cj), j G Z, n-quasiperiodic if there is a sequence 
tj, i G Z, satisfying tjtj+itj^2tj+3 = 1 and such that 

(^j+n = ) bj^n = bj- , Cj^n = Cj- (7) 

Ij+S fj+l 

for each j G Z. 

Note that a sequence tj, j G Z, must be 4-periodic, and it is defined by three parameters, 
e.g., by a := to/ts, (3 := toA2, and 7 := toAi with a(3'j > 0, and hence to = {o:(3'yY^'^. Thus 
the space QSn of n-quasiperiodic sequences has dimension 3n + 3, and {{aj,bj,Cj), j = 
0, n — 1} X (a, (3, 7) are coordinates on it. 

Now we associate a sequence of vectors V^- G M"^, j G Z, and difference equations 

Vj+^ = aj Vj+3 + bj Vj+2 + CjVj+i-Vj, J G Z, (8) 

to each twisted ra-gon Vj := G MP'^,j G Z, with a monodromy M G 5/7(4, M). This 
gives a correspondence between sequences {aj,bj,Cj), j G Z, and twisted n-gons. 
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Proposition 5.9. There is a one-to-one correspondence between twisted n-gons (defined 
up to projective equivalence) and three-parameter equivalence classes in the space QSn of 
n-quasiperiodic sequences {{aj,bj,Cj), j = 0, ...,n — 1} x (a,/3,7). 

If n is odd, there exists a unique n-periodic sequence {aj,bj,Cj) in each class. 

If n = Ap, then the numbers a,P,'-f are projective invariants of a twisted n-gon. 

If n = Ap + 2, then there is one projective invariant: a^yjfi. 

In other words, for odd n the equivalence classes are "directed along" the parameters 
(a, /3, 7) and one can chose a representative with a = /3 = 7 = lin each class. For n = 4p 
the classes are "directed across" these parameters, and hence the latter are fixed for any 
given class. The case n = 4j9 + 2 is in between: in a sense, two of the (a, /3, 7)-parameters 
and one of the (a, 6, c)-coordinates can serve as coordinates on each equivalence class. 

This proposition can be regarded as an analogue of Proposition 4.1 and Remark 4.4 in 
[m for ci = 2. 

Proof. First, we construct the correspondence, and then consider what happens for dif- 
ferent arithmetics of n. For a given n-gon = 0(/c) we construct a sequence of vertices 

e M^ j e Z, in the following way: choose the lifts 0(0) K), 0(1) ^ 14, 0(2) ^ V2 
arbitrarily, and then determine the vectors V^-, j > 2, and Vj, j < 0, recursively using the 
condition det(V^-, V^+i, l^+s) = 1; which follows from equation ([H]). 

By definition of a twisted ra-gon, we have 0(j -|-n) = Mo 0(j) for each j G Z, where M G 
S'L(4, M). Consequently, for each j G Z there exists a number tj, such that Vj+n = tjM^Vj: 
and the matrix M is independent of j. The equation detiVj-^n, l^+i+n, Vj^2+n, V^+s+n) = 1 
implies that tjtj^itj^2tj+3 = 1 and tj^^ = tj for each j G Z. In other words, the whole 
sequence {tj} is determined by ^1,^2, ^3, and then to = l/^i^2^3- Quasiperiodic conditions ([7]) 
follow from the comparison of the equation 

Vj'+4+n = O-j+nYj+3+n + bj-^-nVj^2+n + Cj+nVj+l+n ~ '^j+n 

with equation ([8]). 

Now we rescale the initial three vectors: Vq i— )■ /cqVo, Vi h-> kiVi, V2 ^ ^^2^2, where 
^o^i^2 7^ 0. A different lift of the three initial vectors corresponds to a different sequence 
Vj = kjVj, where the sequence kj,j G Z, must also be 4-periodic and satisfy kQkik2k^ = 1. 
This rescaling gives the action of (M*)'^ on the space QSn of n-quasiperiodic sequences. By 
construction, the corresponding orbits (or equivalence classes) of sequences are in a bijection 
with twisted n-gons. The group (M*)^ acts as follows: 

^ tj{^j/^j+n)y Oij I— >■ ttji^kj / kj^^) , bj H-> bj(kj/kjj^2)-, Cj 1— )■ Cj(kj / kj^i). 
Now we have 3 cases: 

• n is odd. Then the above (R*) ^-action on to? ^i? ^2 allows one to make them all equal to 
1, which corresponds to the constant sequence {tj = 1, j G Z} and a periodic sequence 
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{{aj,bj,Cj), j G Z}. Indeed, e.g., for n = 4p + 3 one has the system of 3 equations: 
toiko/ks) = 1, ti{ki/ki) = 1, ^2(^2/^5) = 1- Since = ko, k^ = fci, and k^ = l/kokik2, 
we obtain a system of 3 equations on the unknowns ko,ki,k2, which has the unique 
solution. 

• n = 4p + 2. One can check that the (M*)^-action does not change the ratio = 

• n = Ap. The (M*) '^-action does not change the three quantities toAs? ^0/^2, ^o/^i- ^ 



Now we can introduce coordinates on the space QSn of n-quasiper iodic sequences using 
Definition 15.51 and quasiperiodic variables {aj,bj,Cj). 

Proposition 5.10. For any n the variables {xj,yj, zj), < j < n — 1 are independent and 
constant on the equivalence classes in QSn, i-e., they are well-defined local coordinates on 
the space Vn = QSn/ ~- 

Proof. It is straightforward to check that the (M*)^-action defined above is trivial on the 
variables {xj,yj,Zj). For instance, 

^ ^ bj+i _^ &j+i(A;j+i/fcj+3) ^ ^ 

^ ajaj+i aj{kj/kj+3) aj+i{kj+i/kj+4) ^ 

The independence of the new variables on QSn follows from that for the original ones. 
Alternatively, it also follows from their local geometric definition (Proposition 15. 7p . □ 

Remark 5.11. In the (a, 6, c)-coordinates for even n some of the (a, /3, 7)-parameters were 
needed to describe the equivalence classes in QSn- In the (x, y, z)-coordinates those param- 
eters are functions of Xj, yj, zf 
i) for n = 4p + 2, 

2p 2 

j-j- x2jy2jZ2j+i _ 0:7 _ 



ii) for n = 4p, 



-p-r X4jX4j+2y4j+2Z4j+l ^ ^ JT l/4j+1^4j ^ ^ JT UijZjj+S ^ 7 

These identities follow from Definitions 15.51 and 15.81 
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6 Algebraic-geometric integrability of the 3D penta- 
gram map 



In this section we complexify the pentagram map and assume that everything is defined over 
C. 

6.1 Scaling transformations and a Lax function in 3D 

Recall that a discrete Lax equation with a spectral parameter is a representation of a dy- 
namical system in the form 

L,,+i(A) = P,+i,(A)L,,(A)i^7/(A), (9) 

where t stands for the discrete time variable, i refers to the vertex index, and A is a complex 
spectral parameter. 

The pivotal property responsible for algebraic-geometric integrability of all pentagram 
maps considered in this paper is the presence of a scaling invariancc. In the 2D case, this 
means the invariance with respect to the transformations aj ajs, bj — )■ bj/s, where s is 
an arbitrary number. In the 3D case, the pentagram map is invariant with respect to the 
transformations aj — )■ ajS, bj — > bj, Cj — > Cjs. In both cases the invariance follows from the 
explicit formulas of the map. Note that formally one can define other pentagram maps by 
choosing the intersection planes in many different ways. However, only very few of these 
maps possess any scaling invariance. Below we derive a Lax representation from the scaling 
invariance. First we do it for odd n, when (a^, bj,Cj),0 < j < n — 1, are coordinates on the 
space Vn- 

Theorem 6.1. The 3D pentagram map on twisted n-gons with odd n admits a Lax repre- 
sentation with the Lax function Lj{\) given by 



/c,/X 


1/A 
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/o 
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Oj/X 








1/A 
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V-1 








) 




^0 





A 





in the coordinates aj,bj,Cj. Its determinant is det Lj = 1/A^. 

Note that we always consider a polygon and the corresponding Lax function at a par- 
ticular moment of time. Whenever necessary we indicate the moment of time explicitly by 
adding the second index "t" to the Lax function (above Lj := Lj^t), while if there is no 
ambiguity we keep only one index. Before proving this theorem we give the following 



24 



Definition 6.2. The monodromy operators Tq^*, Ti^t, T„_i.t are defined as the following 
ordered products of the corresponding Lax functions: 

Tl^t = LQ^tLn-l,tLn-2,t---Ll^ti 

Ti^t = Li^ri-l,tLi-\-n-2,t---Li-\-l^tLi^f, 

Tn-l,t = Ln-2,tLn~3,t---Lo,tLn-l,t^ 

where the (integer) index t represents the moment of time. 

Proof of theorem. First observe that the Lax equation implies that the corresponding 
monodromy operators satisfy 

T,,m(A) = PM(A)T.,(A)P,-i(A), 

i.e., Tj((A) changes to a similar matrix when t — > t + 1, and hence the eigenvalues of the 
matrices Tj t(A) as functions of A are invariants of the map. Conversely, if some function 
Tj_t(A) has this property, then there must exist a matrix -Pi,t(A) (defined up to a multiplication 
by a scalar function) satisfying the above equation. 

How to define such a monodromy depending on a parameter? The monodromy matrix 
associated with the difference equation 

Vj+4 = ajVj+3 + bjVj+2 + CjVj+i - Vj 

is M = NoNiN2...Nn-i, where 

/O -1\ 

10 c, 

^ 10 bj ■ 

\0 1 aj 

For odd n the variables (a^, bj,Cj),0 < j < n — 1, are well-defined coordinates on the space of 
twisted n-gons. These variables are periodic: for any j we have aj+n = = bj, Cj+n = 

Cj. The vectors Vj are quasi-periodic: \^+„ = MVj, and depend on the lift of the points 
from the projective space. This means that the pentagram map preserves the eigenvalues of 
the matrix M, but not the matrix M itself. 

Lemmas 15.31 and 15.41 imply that the pentagram map is invariant with respect to the 
scaling transformations: Therefore, the pentagram map also preserves 
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the eigenvalues of the monodromy matrix M{s) corresponding to the ri-gons scaled by s. 
Namely, we have 



M{s) = No{s)N,{s)N2{s)...N^_,{s), where N^{s) 



/O -1\ 

1 scj 

10 bj 

\0 1 sajj 



Now one can see that the matrix Nj{s) can be chosen as a Lax function. For technical 
reasons (which will be clear later), we define the Lax matrix as Lj^{\) := (g~^Nj{s)g) / s, 
where g := diag(l, s, 1, s), and A := This gives the required matrix Lj{\). □ 



As we mentioned before, the formulas for the pentagram map are non-local in the vari- 
ables {aj,bj,Cj). As a result, an explicit expression for the matrix -Pi,t(A) becomes non-local 
as well. On the other hand, one can use the variables {xj,yj,Zj) (given by Definition 15.51) 
to describe a Lax representation. Their advantage is that all formulas become local and are 
valid for any n, both even and odd. 

Theorem 6.3. For any n the equations for the 3D pentagram map are equivalent to the Lax 
equation 

4,+i(A) = P,+i,(A)4,(A)Pr/(A), 

where 

( 



( 

A(T,-fl + Zi\ 



Vi-idi 



Zi 
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Zi-l 
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1/ 








V 1 + yi^i + Zi 
and the variables pi,ai,Ti, and 9i stand for 
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Pi 







A(Tj 

A 



1 + yi-i + Zi 



-pi \ 

Pi 

1 + yi-2 





Ti 

Oi 



Xi(l + yi + Zi+i) 

Xi-iyi^i{l + yi-2 + Zi_i) 
XiZi-i{l + + Zi){l + yi + Zij^i) ' 

Xiil + yi-2 - ViZi-l + Zi+i + Zi+iyi-2), 

1 + yi-2 + Zi-i 
Ti{l + yi-i + Zi)' 
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Proof. The proof is a long but straightforward verification. 



□ 



Remark 6.4. The Lax functions L and L in the (oj, 6j, q) and (xj, Zi) variables are related 
to each other as follows: 

Li^t = ai+i [h^+iLi^tK) , where K := diag(l, q, 6^, aj). 

6.2 Properties of the spectral curve 

Recall that the monodromy operators Tj ^(A) satisfy the equation 

Tm+i(A) = P,,(A)r,,(A)Pr/(A). 

It implies that the function of two variables R{X, k) = det (Tj ^(A) — k Id) is independent of i 
and t. Furthermore, -R(A, /c) = is a polynomial relation between A and k: R{X, k) becomes 
a polynomial after a multiplication by a power of A. Its coefficients are integrals of motion 
for the pentagram map. The zero set of -R(A, /c) = is an algebraic curve in C^. A standard 
procedure (of adding the infinite points and normalization with a few blow-ups) makes it 
into a compact Riemann surface, which we call the spectral curve and denote by F. In this 
section we explore some of the properties of the spectral curve and, in particular, find its 
genus. 

Definition 6.5. For an odd n define the spectral function R{X, k) as 

R{\k) ■= det (Ti,i(A) - A; Id), 

i.e., using the Lax function in the (a, 6, c)-coordinates from Theorem 16. II The spectral curve 
F is the normalization of the compactification of the curve -R(A, k) = 0. 

We define the integrals of motion Ij, Jj,Gj, < j < q = \n/2\, as the coefficients of the 
expansion 

R{X,k) = k^-k^ W^G,X'-A + k^ ( Xl'^i-^'"'"" ) - ^ ( ) + -^"'^ = 0- 

\i=o / Vi=o / \i=o / 

When n is even, the sequence {aj,bj,Cj), j G Z, is not n-periodic, and the monodromy 
operator Tj_t(A) cannot be defined. One should use the Lax function I/i,t(A) in the {x,y, z)- 
coordinates from Theorem l6.3l to define the monodromy operator Ti^t{X) and the correspond- 
ing spectral curve. 

Namely, first note that the integral of motion Jq has the following explicit expression: 

n-l /n-1 \ ^1/4 

/o = n = In ) 

i=0 \j=0 / 
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Definition 6.6. For any n (either even or odd), the spectral function is 

R{\, k) = R{X, kIo)/Io, where ^(A, k) := det (Ti,t(A) - A;Id) , 

and the monodromy operator Tj((A) is defined using the Lax function Li^t{X) from Theo- 
rem [631 

The spectral function R{X, k) defined this way coincides with det (Tj ((A) — A; Id) for odd 
ra, since Tj ^ = Iq [hl^Ti^thi) , see Remark 16.41 It is convenient to have such a unified definition 
for computations of integrals of motion. 

Theorem 6.7. For generic values of the integrals of motion Ij, Jj,Gj, the genus g of the 
spectral curve T is g = 3q for odd n and g = 3q — 3 for even n, where q = \n/2\. 

To prove it, we first describe the singularities of -R(A, /c) = by considering the formal 
series solutions (the so-called Puiseux series). 

Lemma 6.8. If n is even, the equation R{X, k) = has 4 formal series solutions at X = 0: 



Oi : k, = ^-^X + 0{X'' 
-'0 



^2,3 '■ ^2,3 ^ + ^ ( A^ j ' '"^^'^^'^ ^* satisfies G^kl - J^k^ + Jq = 0, 



and 4 solutions at X = oo 



W^i,2,3,4 : A;i,2,3,4 = ^ + ^ ()^) ' ^here k^ - G,kl + Jgkl - I,k^ + 1 = 0. 



If n is odd, the equation R{X, k) = has 4 formal series solutions at X = 0: 

Oi : fci = i - ^A + OiX^' 
-'0 



A"/2 2GoA("-i)/2 VA("-2)/2 ; ' 

Gi C 
— H 

A" ' A""i A*^ 



and 4 solutions at X = oc: 

The remaining coefficients of the series are determined uniquely. 
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Proof of lemma. One finds the series coefficients recursively by substituting the series 
into the equation R{X, k) = 0, which determines the spectral curve. □ 



Now we can complete the proof of Theorem 16.71 

Proof of theorem. As follows from the definition of the spectral curve F, it is a ramified 
4- fold cover of CP^, since the 4 x 4- matrix Tj^t (or Ti^t) has 4 eigenvalues. By the Riemann- 
Hurwitz formula the Euler characteristic of F is x(r) = 4x(CP ) — u = 8 — u, where u is the 
number of branch points. On the other hand, x(r) = 2 — 2g, and once we know u it allows 
us to find the genus of the spectral curve F from the formula 2 — 2g = 8 — u. 

The number u of branch points of F on the A-plane equals the number of zeroes of the 
function dkR{X, k) aside from the singular points. The function dkR{X, k) is meromorphic 
on F, therefore the number of its zeroes equals the number of its poles. One can see that for 
any n the function dkR{X, k) has poles of total order 9n at 2; = 0, and it has zeroes of total 
order 6n at z = 00. Indeed, substitute the local series from Lemma 16.81 to the expression 
for dkR{X,k). (E.g., for n = 2q ai one has k ~ A~". The leading terms of dkR{X,k) 
for the pole at A = are k^, A;^A~", kX~'^~^, A~^". The first two terms, being of order A~^", 
dominate and give the pole order of 3n = 6q.) The corresponding orders of the poles and 
zeroes of dkR{X, k) on F are summarized as follows: 
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For instance, for n = 2q this gives the total order of poles: Aq + Aq + Aq + Qq = 18q = 9n, 
while the total order of zeroes is 4 x 3g = 12q = 6n. 

Therefore, the number of zeroes of dkR{X, k) at nonsingular points A 7^ {0, 00} is u = 
9n — 6n = 3n, and so is the total number of branch points of F in the finite part of the (A, k) 
plane. The difference between odd and even values of n arises because F has 2 additional 
branch points at A = 00, and 1 branch point at A = for odd n, i.e., u = 3n + 3. 

Consequently, one has 2 — 2g = 8 — u with u = 3n for even n and u = 3n + 3 for odd 
n. The required expression for the genus g follows: g = 3q — 3 for n = 2q and g = 3q for 
n = 2q + l. □ 

Remark 6.9. Now we describe a few integrals of motion using the coordinates {aj, bj, cj), < 
j < n — 1, when n is odd. The description is similar to that in the 2D case (cf. Section 5.2 
and Proposition 5.3 in [llj). Consider a code which is an ordered sequence of digits from 1 
to 4. The number of digits in a code is p, q, r, t, respectively. The code is called "admissible" 
if p + 2q + 3r + At = n. The number p + r is called its "weight." Each code expands in a 
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"word" of n characters in the following way: 1,2,3,4 are replaced by "a" , "*b" , "**c" , ^ 
respectively. Now we label the vertices of an ra-gon by 0, 1, ...,n — 1, and associate letters in 
a word to them keeping the order. We obtain one monomial by taking the product of the 
variables aj,6j, or q that occur at the vertex i. The letter "*" corresponds to "1". The sign 
of the monomial is (—1)*. Next step is to permute the numbering of the vertices cyclically 
and take the sum of the monomials. Note, however, that if, for example, n = 9, then the 
code "333" corresponds to the sum cqCsCg + C1C4C7 + C2C5C8 without the coefficient 3. Finally, 
we repeat this procedure for all admissible codes of weight p + r and denote the total sum 
by Ip+r- Additionally, we define the sum Gp+r by substituting Oj — )■ q — )■ aj_i in Ip+r- 

Consider, for example, the case n = 7. Then all admissible codes of weight 1 are 
142, 124, 1222, 34, 223. The corresponding sum is 

Ji = ^ (-ai6o - 0560 + a56o&2&4 - cq + co62&4)- 

cyclic 

Proposition 6.10. For odd n one has 

III h = In-2i 5 

n-1 

Gq = Y\ Cj ) Gq = Gi, Gi = Gn-2i ■ 
j=0 

Proof. The proof is analogous to the proof of Proposition 5.3 in [TT]. □ 



n 

1=0 



6.3 The spectral curve and invariant tori 

The spectral curve is a crucial component of algebraic-geometric integrability. Below we 
always assume it to be generic. (As everywhere in this paper, "generic" means values of 
parameters from some Zariski open subset in the space of parameters.) It has a natural 
torus, its Jacobian, associated with it. It turns out that one can recover a Lax function 
from the spectral curve and a point on the Jacobian, and vice versa: in our situation this 
correspondence is locally one-to-one. The dynamics of the pentagram map becomes very 
simple on the Jacobian. In this section, we construct this correspondence and describe the 
dynamics of the pentagram map. 

Definition 6.11. A Floquet-Bloch solution ipi^t of a difference equation ipi+i^t = Li,t'ipi,t is an 
eigenvector of the monodromy operator: 
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The normalization ^^=1 V'o,o,j = 1 (i-G-? the sum of all components of the vector tpo^ is equal 
to 1) determines all vectors ipi^t with i,t > uniquely. Denote the normalized vectors ipi^t by 

ipi^t, i-e., ipi^t = i^i^t/ (j2'j=ii'i,t,j^ ■ (The vectors ipo^ and -i/^o.o are identical in this notation.) 

We also denote by Di t the pole divisor of tpi^t on T. 

Remark 6.12. We use the Lax function Lj ^ and the monodromy operator Tj j in the above 
definition to allow for both even and odd n. In the case of odd n one can instead employ 
the Lax function Lj j and the monodromy operator Tj j, while all the statements and proofs 
below remain valid. 

Theorem 6.13. A Floquet-Bloch solution ipi^t is a meromorphic vector function on T. 
Generically its pole divisor Di^t has degree g + 3. 

Proof. The proof of the fact that the function ipi^t is meromorphic on the spectral curve F, 
as well as that its number of poles is deg Di^t = ^/2, is identical to the proof of Proposition 
3.4 in [15]. The number u of the branch points of F is different: in Theorem 16.71 we found 
that 2 — 2g = 8 — u, where g is the genus of the spectral curve. This implies the required 
expression: deg Dit = + 3. □ 

Definition 6.14. Let J(F) be the Jacobian of the spectral curve F, and [-Do,o] is the equiv- 
alence class of the divisor -Do.o, the pole divisor of ipo^, under the Abel map. The pair 
consisting of the spectral curve F (with marked points Oi and Wi) and a point [-Do,o] £ J(X) 
is called the spectral data. The spectral map S associates to a given generic twisted n-gon in 
CP^ its spectral data (F, [Do,o]). 

The algebraic-geometric integrability is based on the following theorem. 

Theorem 6.15. For any n, the spectral map S defines a bijection between a Zariski open 
subset of the space Vn = {{xi,yi, Zi), < i < n — 1} and a Zariski open subset of the spectral 
data. 

Corollary 6.16. For odd n, the spectral map S defines a bijection between a Zariski open 
subset of the space Vn — C^" = {(cti, q), < i < n — 1} and a Zariski open subset of the 
spectral data. 

Proof of Corollary 16.161 The statement follows from Theorem 16.151 and Definition 15.51 
relating the coordinates (a;^, Zi) and (aj, 6j, q) for odd n. □ 

The proof of Theorem l6.15l is based on Proposition 16 . 1 71 (which completes the construction 
of the direct spectral map) and Proposition 16.181 (an independent construction of the inverse 
spectral map), which are also used below to describe the corresponding pentagram dynamics. 
It will be convenient to introduce the following notation for divisors: Opq := Op + Oq and 
Wpq ■= Wp + Wq (e.g., Oi2 := Ol + O2). 
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Proposition 6.17. The divisors of the coordinate functions ipi.t,!, V'i.t,4 for < i < n — 1 
and any integer t satisfy the following inequalities, provided that their divisors remain non- 
special up to time t: 
For odd n one has 

• (^i,t,i) >-D + 02- z023 + + l)Wi2 - t{Wi2 - Ols); 

• (V^M,2) > + (1 - i)02z + iWi2 - t{Wi2 - Ols); 

• (^M,3) >-D- 1O2Z + + l)Wi2 - t{Wi2 - O13); 

• {'ipiM) >-D + 02 + {l- OO23 + iWi2 - t{Wi2 - O13); 
For even n one has 

• (^.,,1) >-D + 02+ f'^^^^ J W,2 + L^^^J W^34 - L^J O24 - L^J O34 + tOu; 

• {i'^,t,2) >-D+ L^^^J + L^J Wu - 02A - L^J O34 + to,,- 

• (^M,3) > + L^^^J W^12 + L^^^J W^34 - L^J O34 - [\\ O24 + tOu; 

• (V'.m) > -/^ + O2 + L^^^J w^i2 + L^J w^34 - L^J O34 - L^J O24 + to,,- 

where D = Do,o corresponds to the divisor att = and is an effective divisor of degree g + 3, 
while [x\ is the floor (i.e., the greatest integer) function of x. 

Proof. The proof is a routine comparison of power expansions in A at tlie points Op, Wq 
for ki and Li^t and is very similar to the proof of Proposition 3.10 in the 2D case in [15], 
although the 3D explicit expressions are more involved. See more details in Appendix 19.31 

□ 



Proposition 6.18. For any n, given a generic spectral curve with marked points and a 
generic divisor D of degree (7 + 3 one can recover a sequence of matrices 



/ -i\ 

Xxiyi 1 



-1 





\ 







Xxi 1 J 



for < i < n — 1 and any t. 

We describe the reconstruction procedure and prove this proposition in Appendix 19.31 
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Proof of Theorem 16.151 The proof consists of constructions of the spectral map S and 
its inverse. The spectral map was described in Definition 16.141 based on Theorem 16.131 We 
comment on an independent construction of the inverse spectral map now. 

Pick an arbitrary divisor D of degree g+3 in the equivalence class [-Do,o] ^ <^(r) and apply 
Proposition l6.18[ "A Zariski open subset of the spectral data" is defined by spectral functions 
which may be singular only at the points Oi, Wi and by such divisors [D] = [-Do,o] ^ J(X) 
that all divisors in Proposition 6.17 with 0<i<n — lupto time t are non-special. □ 

The next theorem describes the time evolution of the pentagram map in the Jacobian of P. 
The difference between even and odd n is very similar to the 2-dimensional case. Combined 
with Theorem I6.15[ it proves the algebraic-geometric integrability of the 3D pentagram 
map. (It also implies that it is possible to obtain explicit formulas of the coordinates of the 
pentagram map as functions of time using the Riemann ^-functions.) 

Theorem 6.19. The equivalence class [-Dj,t] € J(P) of the pole divisor Dn of ipi^t has the 
following time evolution: 

• when n is odd, 

[Di,t] = [Do,o - tOn + i023 + it- t)Wi2l 



when n is even, 



where deg Di^t = fl' + 3 and \_x\ is the floor function of x, and provided that spectral data 
remains generic up to time t. 

For an odd n this discrete time evolution in J(P) takes place along a straight line, whereas 
for an even n the evolution goes along a "staircase" (i.e., its square goes along a straight 
line). 

Proof. The vector functions ipi^t with i,t are not normalized. The normalized vectors 
are equal to tpi^t = 'ipi,t/ fi,t, where fi^t = J2'j=i'^iAj- Proposition 16. 171 implies that the divisor 
of each function fn is: 



for odd n, 



for even n, 



12, 
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Since the divisor of any meromorphic function is equivalent to zero, the result of the theorem 
follows. The staircase dynamics is related to alternating jumps in the terms — t + 1)/2J 
and [(z — t)/2\ as t increases over integers. □ 

Note that although the pentagram map preserves the spectral curve, it exchanges the 
marked points. The "staircase" dynamics on the Jacobian appears after the identification 
of curves with different marking. One cannot observe this dynamics in the space of twisted 
polygons Vn, before the application of the spectral map. 

7 Ramifications: closed polygons and symplectic leaves 
7.1 Closed polygons 

Closed polygons in CP^ correspond to the monodromies M = ±ld in SL(A,C). They form 
a subspace of codimension 15 = dimS'L(4,C) in the space of all twisted polygons Vn- The 
pentagram map on closed polygons in 3D is defined for n > 7. 

Theorem 7.1. Closed polygons in CP^ are singled out by the condition that either {\,k) = 
(1,1) or {X,k) = (1,-1) is a quadruple point of T . Both conditions are equivalent to 9 
independent linear constraints on Ij, Jj, Gj. Generically, the genus ofT drops to g = 3q — 9 
when n is even, and to g = 3q — 6 when n is odd, where q = [n/2\. The dimension of 
the Jacobian J{T) drops by 6 for closed polygons for any n. Theorem \6.13\ holds with this 
genus adjustment, and Theorems \6.15\ and \6.1^h hold verbatim for closed polygons (i.e., on 
the subspace of closed polygons Cn C Vn)- 

Proof. For a twisted ra-gon its monodromy matrix at a moment t is equal to To^t(l) in 
the (a, 6, c)-coordinates or to To^((l) in the (x, ?/, 2;)-coordinates. An n-gon is closed if and 
only if To,t(l) = Id or To,t(l) = -Id (respectively, To,t(l) = /qM or To,i(l) = -/old). 
For our definition of the spectral function, either of these conditions, To^i(l) = ±Id or 
T^ tiX) = ±Iq Id, implies that (A, A;) = (1, ±1) is a self-intersection point for F. 
The algebraic conditions implying that (1, ±1) is a quadruple point are: 

• R{l,±l) = 0, 

• M(l,±l) = M(1,±1) = 0, 

. dlRil, ±1) = dlRil, ±1) = dlR{l, ±1) = 0, 

. 9|i?(l, ±1) = dlRil, ±1) = dl.Ril, ±1) = dl,R{l, ±1) = 0. 

However, the function R{X, k) is special at the points (1, ±1), because the following relation 
holds: 

i?(l, ±1) = ±dkRil, ±1) - ^dlRil, ±1) ± ^dlRil, ±1). 
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Consequently, the above 10 conditions are equivalent to only 9 independent linear equations 
on Ij, Jj, Gj, < j < q. 

The proofs of Theorems 16.151 and 16.191 apply, mutatis mutandis, to the periodic case. 
To define the Zariski open set of spectral data for closed polygons, we confine to spectral 
functions that can be singular only at the point {X,k) = (1,1) or (1,-1) in addition to 
singularities at Oi and Wi and use the same restrictions on divisors D as in the proof of 
Theorem 16.151 

In the periodic case we also have to adjust the count of the number u of branch points 
of r and the corresponding calculation for the genus g of F, cf. Theorem I6.13[ Namely, as 
before, the function dkR{X, k) has poles of total order 9n over A = 0, and zeroes of total 
order 6n over A = oo. Now since R{X,k) has a quadruple point (1,±1), dkR{X,k) has a 
triple zero at (1,±1). But A = 1 is not a branch point of F. Consequently, dkR{X,k) has 
triple zeroes on 4 sheets of F over A = 1. The Riemann-Hurwitz formula is 2 — 2(7 = 8 — i/, 
where the number of branch points for even n is = 9n — 6n — 12 = 3n — 12, while for odd n 
it is = 9n — 6n — 12 + 3 = 3n — 9. Therefore, we have g = 3q — 9 for even n, and g = 3q — 6 
for odd n. □ 



7.2 Invariant symplectic structure and symplectic leaves 

It was proved in [15j that in the 2D case an invariant symplectic structure on the space of 
twisted polygons P„ provided by Krichever-Phong's universal formula [5l E] coincides with 
the inverse of the invariant Poisson structure found in [llj when restricted to the symplectic 
leaves. We show that in 3D the same formula also provides an invariant symplectic struc- 
ture defined on leaves described below. While we do not compute the symplectic structure 
explicitly in the coordinates (aj,5j,Cj) or {xi,yi,Zi) due to complexity of the formulas, the 
proofs are universal and applicable in the higher- dimensional case of CP^ as well. Finding 
an explicit expression of the symplectic structure or of the corresponding Poisson structure 
is still an open problem. 

Definition 7.2 (OE]). Krichever-Phong's universal formula defines a pre- symplectic form 
on the space of Lax operators, i.e., on the space Vn- It is given by the expression: 



The matrix \E'o := \E'o,t(A) is composed of the eigenvectors ipo,t on different sheets of F over 
the A-plane, and it diagonalizes the monodromy matrix Tq := To,f(A). (In this definition we 
drop the index t, because all variables correspond to the same moment of time.) 

The leaves of the 2-form u are defined as submanifolds of Vn, where the expression 
6 Ink (dX/X) is holomorphic. The latter expression is considered as a 1-form on the spectral 
curve F. 




A=0,oo 
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Proposition 7.3. For even n the leaves are singled out by 6 conditions: 



5Io = 5Iq = 5Go = 5Gq = 5Jo = 5Jq = 0; 
For odd n the leaves are singled out by 3 conditions: 

SGq = SIq = 5Jq = 0. 

Proof. These conditions follow immediately from the definition of the leaves and Lemma [6l8l 
For example, at the point Oi we have 

This 1-form is holomorphic in A if and only if SIq = 0. Similarly, we obtain 6{Iq/Gq) = 
at the point O2 for odd n. (One has to keep in mind that the local parameter around this 
point is A^/^) □ 

Remark 7.4. The definition of a presymplectic structure u on Vn uses \fo cind Tq and 
hence relies on the normalization of \E'o. When restricted to the leaves from Proposition I7.3[ 
the 2-form u becomes independent of the normalization of the Floquet-Bloch solutions. 
Additionally, the form u becomes non-degenerate, i.e., symplectic, when restricted to these 
leaves, as we prove below. The symplectic form is invariant with respect to the evolution 
given by the Lax equation. The proof is very similar to that of Corollary 4.2 in |8] (cf. [5| 
for other proofs). 



Theorem 7.5. The rank of the invariant 2-form u restricted to the leaves of Proposition \L3 
is equal to 2g. 

Proof. Since the 1-form 5hik dX/X is holomorphic on F, it can be represented as a sum of 
the basis holomorphic differentials: 

d\ ^ 

Sink — = Y^SUidu„ (10) 

i=l 

where g is the genus of F. The coefficients Ui can be found by integrating the last expression 
over the basis cycles of Hi{T): 

Ui = (f \nk^ . 

According to formula (5.7) in [7], we have: 

9+3 

tu = 6 In k{pi) A 6 In A(pi), 

i=l 
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where the points pi E T, 1 < i < g + 3, constitute the pole divisor Dq.o of the normahzed 
Floquet-Bloch solution ^0,0 • 

After rearranging the terms, we obtain: 

'9+3 



UJ 



where 



9+\^ CPs 

s=l 

are coordinates on the Jacobian J(r). The variables Ui and ipi are natural Darboux coor- 
dinates for u, which also turn out to be action-angle coordinates for the pentagram map. 
(The latter follows from the general properties of the Krichever-Phong universal form for a 
given Lax representation, cf. [3 |6].) 

Let us show that the functions Ui are independent. Assume the contrary, then there 
exists a vector v on the space Vn, such that 6Ui{v) = for all i. Then it follows from ( !T0|) 
that dyk = 0. After applying the operator to R{X,k), we conclude that k satisfies an 
algebraic equation of degree 3, which is impossible, since F is a 4-fold cover of the A-plane. 
□ 



Remark 7.6. In more details, there are the following two cases: 

• even n = 2q. The dimension of the space Vn is 6g. The codimension of the leaves is 
6. Therefore, the dimension of the leaves matches the doubled dimension of the tori: 
2g = 6q — 6. 

• odd n = 2g + 1. The dimension of the space Vn is Qq + 3. The codimension of the 
leaves is 3. Again, the dimension of the leaves matches the doubled dimension of the 
tori: 2g = 6q. 

The algebraic-geometric integrability in the complex case implies Arnold-Liouville inte- 
grability in the real one. Indeed, the pre-symplectic form depends on entries of the mon- 
odromy matrix in a rational way, since it is independent of the permutation of sheets of the 
spectral curve F. Therefore, its restriction to the space of the real n-gons provides a real 
pre-symplectic structure. One obtains invariant Poisson brackets on the space of polygons 
Vn by inverting the real symplectic structure on the leaves, while employing invariants of 
Proposition 17.31 as the corresponding Casimirs. 

Problem 7.7. Find an explicit formula for an invariant Poisson structure with the above 
symplectic leaves. 
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8 A Lax representation in higher dimensions 



The origin of the integrabihty of the pentagram map is the presence of its scahng invariance. 
Assume that gcd{n,d+ 1) = 1. The difference equation ([T]) 

Vj+d+i = djA^j+d + aj^d-iVj+d-i + ••• + CLj^y-j+i + {-^YVj 

allows one to introduce coordinates aj^i,aj^2, ■■■o,j,d, < j < n — 1, on the space of twisted 
n-gons in any dimension d. 

Proposition-conjecture 8.1. (The scaling invariance) The pentagram map on twisted 
n-gons in CW^ is invariant with respect to the following scaling transformations: 

• for odd d = 2k + 1 the transformations are 

o-j,! sttj^i, Oj^s — )■ sOj^s, 5 — )■ sa^ s, ... ,aj^d saj^d, 
while other coefficients aj^2i with I = 1, x do not change; 

• for even d = 2x the transformations are 

Oij 2 — ^ S(lj 2^ ^jA — ^ ^ '^ii4' ■■■ ' — ^ ^jA 

for all seC*% 

Proof. In any dimension d the pentagram map is a composition of involutions a and /3, 
see Section 15.11 (More precisely, a is not an involution for even d, but its square is a shift 
in the vertex index, see [H] for the 2D case.) One can prove that the involution a : Vi Wi 
in any dimension has the form 

Wj+d+i = (-1)^^+^(0^,11^^-+^ + a^^2Wj+d--i + ... + a^,dWj+i - Wj) , 

where * stands for the first index, which is irrelevant for the scaling (Lemma 15.31 proves the 
case d = 3). 

We call this Proposition-conjecture because the proof of an analog of Lemma [5.41 (for the 
map P) in higher dimensions is computer assisted. One verifies that for a given dimension d 
the coefficients consist of the terms that are consistent with the scaling. □ 

We obtained explicit formulas, and hence a direct (theoretical) proof of the scaling in- 
variance for the pentagram maps up to dimension d < 6. This bound is related to computing 
powers to produce explicit formulas and might be extended. However, we have no general 
purely theoretical proof valid for all d and it would be very interesting to find it. 

"'Wc thank G.Mari-Bcffa for correcting an error in the scahng for even d in the first version of this 
manuscript, as weU as in the short version |16) . This error related to numerics with a different choice of 
vertices for the diagonal planes leads to another system, different from Tp^r, which also turns out to be 
integrable and will be discussed elsewhere. 
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Problem 8.2. Find a general proof of the scaling invariance of the pentagram map in any 
dimension d. 



Theorem 8.3. The scale-invariant pentagram map on twisted n-gons in any dimension d is 
a completely integrable system. It is described by the Lax matrix 



/ 







\ 



aj-2 



where D{X) is the following diagonal matrix of size d x d: 

• for odd d = 2x + 1, one has D{X) = diag{X, 1, A, 1, A); 

• for even d = 2h, one has one has -D(A) = diag{l, A, 1, A, . 



1,A). 



Proof sketch. By using the scahng invariance of the pentagram map, one derives the 
Lax matrix exactly in the same way as in 3D, see Section 16.11 Namely, first construct the 
(d+l) X ((i+l)-matrix Nj{s) depending on our scaling parameter s, and then use the formula 
L~^(A) = {g~'^Nj{s)g) /s™ with a suitable choice of the diagonal {d + 1) x {d + l)-matrix g 
and an appropriate function of the parameter s. 

For odd d = 2x+ 1, we have g = diag(l, s, 1, s, 1, s), m = 1, and A = whereas for 
even d = 2x, we have g = diag{l, s''^, s, s^''^, s"^, s'^"^, s"^, s'^), m = x, and A = s~'^~^. 
The Lax representation with a spectral parameter is constructed as we described above. 

Using the genericity assumptions similar to those used in the 2D and 3D cases, one 
constructs the spectral map and its inverse, which is equivalent to algebraic-geometric inte- 
grability of the pentagram map. Coefficients of the spectral curve form a maximal family 
of first integrals. Along with a (pre)symplectic structure defined by the Krichever-Phong 
formula, this provides the Arnold-Liouville integrability of the system on the corresponding 
symplectic leaves in the real case. □ 



The scaling parameter has a clear meaning in the continuous limit: 

Proposition 8.4. For any dimension d the continuous limit of the scaling transformations 
corresponds to the spectral shift L L + X of the differential operator L. 

Proof. In 2D this was proved in [TTj. A continuous analog of the difference equation ([1]) is 

G{x + ((i + l)e) = ad(x, e)G{x + de) + ... + ai(x, e)G{x + t) + {-ifGix), 

where G{x) satisfies the differential equation (|2]) with a differential operator L of the form (|5]). 
Using the Taylor expansion for G{x+ie) and the expansion aj{x, e) = a°(a;) + eaj(a;) + we 
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obtain expressions of a], in terms of the coefficients of L, i.e., in terms of functions Uj{x) and 
their derivatives. We find that the terms a° are constant, = for all k, while a], for i > 2 
are linear in Ud-i+i and differential polynomials in the preceding coefficients Ud-i, Ud-i+2- 
The scaling parameter also has an expansion in e: s = ro+eri+e^r2+.... We apply it to the 
coefficients ai{x,e) and impose the condition that and = are fixed, similarly to [Ilj. 
By term-wise calculations (different in the cases of even and odd d and using the "triangular 
form" of the expressions for al), one successively obtains that tq = 1, ti = ... = = 0, i.e., 
s must have the form s = 1 + Td+ie'^^^ + 0{e'^~^'^). Its action shifts only the last term of L: 
Uq ^ Uq + const ■ Trf+ie'^^^, i.e., it is equivalent to the spectral shift L — )■ L + A. □ 

Note that the spectral shift commutes with the KdV flows. Indeed, d/dt{L + A) = 
d/dt L = [Q2, L] = [Q2, L + X], since Q2{L) := d'^ + = Q2{L + A) for operators L of 

degree d + 1 > 3. Equivalently, the pentagram map commutes with the scaling transforma- 
tions in the continuous limit. 

9 Appendices 

9.1 Continuous limit in the 3D case 

In this section we present explicit formulas manifesting Theorem |X] on the continuous limit 
of the 3D pentagram map. Consider a curve G{x) in given by the differential equation 

G"" + u{x)G" + v{x)G' + w{x)G = 

with periodic coefficients u{x),v{x),w{x). To find the continuous limit, we fix e and consider 
a plane Pe{x) passing through the three points G{x — e), G{x), G{x + e) on this curve. We 
are looking for an equation of the envelope curve L^{x) for these planes. 
This envelope curve L^{x) satisfies the following system of equations: 

det\G{x),G{x + e),G{x-e),L,{x)\ = 
det 1^(0;), G{x + e), G{x - e), L'^{x) \ = 
det|G(x),G(x + e),G(x-e),L';(x)| =0. 

By considering the Taylor expansion and using the normalizations det\L^, L'^, L", L"'\ = 1 
and det |G, G', G", G"'\ = 1 we find that 

L^Xx) = G{x) + ^ (g'\x) + ^G{x)) + 0{e') (11) 
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as e — )■ 0. Now, the equation L'^" + + v^L'^ + w^L^ = imphes that: 

v, = v + ^i2w' + v" - uu' - 2m'") + C(e^), 
6 

w, = w + ^{2w" - uu" - vu! - u!'") + C(e^). 
These equations describe the (2, 4)-equation in the Gelfand-Dickey hierarchy: 

{u = 2v' - 2u'\ 
v = v" + 2w' - uu' - 2u!\ 
111 
// / I! Hit 
w = w vu uu u , 
2 2 2 ' 

where L = d'^ + ud"^ + vd + w and Q2 = (-^^^^^)+ = d'^ + 

Remark 9.1. A different choice of the points defining the plane Pe{x) on the original curve 
leads to the same continuous limit. For instance, the choice of G{x — 3e), G{x + e), G{x + 2e) 
results in the same expression for L^{x), where in ffTTj) instead of the coefficient one has 
7e^/6. This leads to the same evolution of the curve G with a different time parameterization, 
cf. Remark 14.41 



9.2 Higher terms of the continuous Umit 

Recall that in the continuous limit for the pentagram map in MF'^ the envelope for osculating 
planes moves according to the {2,d + l)-KdV equation (Theorem 14.51) . This evolution is 
defined by the e^-term of the expansion of the function L^{x). 

The same proof works in the following more general setting. Let L be a differential 
operator of order d + 1 and G a non-degenerate curve defined by its solutions: LG = 0. 

Proposition 9.2. Assume that the curve G evolves according to the law G = QmG, where 
Qm ■= (L™'/*^^^^))4. is the differential part of the mth power of the operator Q = L^/(^+'^\ 
Then this evolution defines the equation L = [Qm,L], which is the {m,d + l)- equation in the 
corresponding KdV hierarchy of L. 

Furthermore, one can define the simultaneous evolution of all terms in the e-expansion 
of L^{x) using the following construction. For the pseudodifferential operator Q := L^/^'^+^'i 
consider the formal series exp(eQ) := 1 + eQ + yQ^ + ••• and take its differential part: 

(exp(eg))^ = \^l + eQ + -Q^ + ... j = 1 + eQi + -Q2 + - = —\^^ ■ 
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For each power of e this is a multiple of the differential operator Qm, which is the differential 
part of the mth power of the operator Q = 1^1^'^+^). 

Corollary 9.3. The formal evolution equation G = {exp{eQ))^G corresponds to the full 
KdV hierarchy L = [(exp(eQ))_^ , L], where the operator L is of order d+1 and the {m,d+l)- 
equation corresponds to the power e™". 

A natural question is which equations of this hierarchy actually appear as the evolution 
of the envelope L^{x). Recall that only even powers of e arise in the expansion of the function 
L^[x) for the continuous limit of the pentagram map. The e^-term gives the (2, d + l)-KdV 
equation. It turns out that the e^-term in the continuous limit of the 2D pentagram map 
results in the equation very similar to the (4, 3)-equation in the KdV hierarchy (which is a 
higher-order Boussinesq equation). Although the numerical coefficients in these differential 
equations are different, one may hope to obtain the exact equations of the KdV hierarchy 
for different m by using an appropriate rescaling. This allows one to formulate 

Problem 9.4. Do higher (m, d + l)-KdV flows appear as the e'^-terms in the expansion of 
the envelope L^{x) for the continuous limit of the pentagram map for any even m > 2? 



9.3 Bijection of the spectral map 

In this appendix we sketch the proof of Proposition 16.171 and prove Proposition 16.181 which 
allows one to reconstruct the L-matrix from spectral data, and hence complete the proof of 
Theorem 16.151 on the spectral map. 

Proposition 9.5. (= Proposition 16.181) For any n, given a generic spectral curve with 
marked points and a generic divisor D of degree (7 + 3 one can recover a sequence of matrices 



( ° 








-1\ 










1 





Zi 





1 


I 





Xxi 


1/ 



for < i < n — 1 and any t. 



Proof. Without loss of generality we describe the procedure to reconstruct the matrices 
Li{\) := Z,,o(A), for < i < n - 1 and t = 0. 

1. First, we pick functions ipij '■= V'i.Oj for < i < n — 1, 1 < j < 4, and t = 0, satisfying 
Proposition 16.171 Note that according to the Riemann-Roch theorem, the functions 
ipi^i and ipi^4 are deflned up to a multiplication by constants, whereas the functions ipi^2 
and ipi^s belong to 2-dimensional subspaces. The functions ipi^i and ipi^4 belong to the 
same subspaces. We pick the pairs of functions ipi^i, ipi^s and 'ipi^2, ipi,A to be linearly 
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independent. Observe that any sets of functions ipiA satisfying Proposition 16.171 

are related by gauge transformations ipi — )■ g~'^ipi, where 



/AO o\ 

Bi E, 

Fi Ci 

\0 OA/ 



9i 9i+ny 



and ipi stands for ipi = {ipi,iy ■ We also define ipn to be ipn = lokipo for any n 

in {x, y, z)-variables. 

2. We find the unique matrix L'^ satisfying the equation ipi = (L^'^V'j+i- 



/ ' 

t,,3 



3. One can check that there exists the unique choice of the matrices Qi, < i < n — 1, 
such that the equality Lj(A) = gi+iL'-(X)g~^ is possible. The latter is equivalent to the 
following system of equations (0 < i < n — 1): 



D 



i+1 



D 



1; 



i+1 



A+i A+i 

These equations decouple and may be solved explicitly. One only needs to check the 
solvability of n equations — — = 1, 0<i<n — 1, for n variables A, < i < n — 1. 



0. 



i+1 



A non-trivial solution exists provided that 11^=0^ '''«:4 ~ depends on an arbitrary 
constant, which corresponds to multiplication of all matrices gi by the same number 
and does not affect the Lax matrices. One can check that 



f-iA 



V^^+l,4(Ql) 
i^iA{0,) 



n-1 



and JJ^ 



'iA 



1=0 



^n,4(Ql 



By using Lemma [6.81 we find the value k{Oi) = I/Iq as required. Now the remaining 
variables Ai, Bi,Ci, Ei, Fi, < i < n — 1, are uniquely determined. 



□ 
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Corollary 9.6. For odd n, given a generic spectral curve with marked points and a generic 
divisor D one can recover a sequence of matrices 



with Q < i < n — 1 and any t. 



/O -1\ 



A 



-1 



1 
\0 A aj / 



Proof. The statement follows from Proposition 16.181 and the fact that {aj,bj,Cj), <i < 
n — 1, are coordinates on the space Vn for odd n. □ 

We complete the exposition with a sketch of the proof for Proposition 16.171 for even n 
(the case of odd n is similar). 

Proposition 9.7. (= Proposition 16.171 ) For even n, the divisors of the coordinate func- 
tions V'j,*,!, '?/'i,t,4 for Q < i < n — 1 and any integer t satisfy the following inequalities, 
provided that their divisors remain non-special up to time t: 

• (Am) >-D + 0,+ C-^^\ W^i2 + C~l~^^ \ Wu - L^J - L^J O34 + tOu; 

• {i^^,t,2) >-D+ [ '~\^^ \ w,2 + L^J w^34 - L^J O24 - L^J O34 + to,,- 

• (^M,3) > + L^^^J + L^^^J ^34 - L^J O34 - L^J O24 + tOu; 

• (V'.m) >-d + 02+ L^^^J + L^J w^34 - L^J O34 - L^J O24 + tOi4; 

where D is an effective divisor of degree g + and \_x\ is the floor function of x. 



Proof. First, we prove these inequalities for t = and < i < n — 1. For illustration 
we find the multiplicities of the components of the vector ipi^ at the point O2, while other 
points can be treated in a similar fashion. We employ the matrices Lj ^ in the coordinates 

•^i) Ui) ^i- 

Notice that a cyclic permutation of indices (n — 1, n — 2, 1, 0) changes the monodromies 
Tj — )■ Tj+i and the Floquet-Bloch solutions Tpi — )■ V^j+i. For even n, it also permutes ipi{02) ^ 
ipii^Os) and Wu f-)- VF34, i.e., the corresponding pairs of the vectors ipi at the points {Wi, W2) 
and (VI/3, W4) are swapped. 
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Using the asymptotic expansion of To,t(A) at A = 0, the definition of the Floquet- 
Bloch solution, and the normahzation condition, one can show that t/^q.o = {0{)^), 0(A), 1 + 
0(A), O(A))"^ as A — )■ at the point O2. Since 



(\ 10 



y\ yi 

\o 








as A — 7- 0, 



and 'i/'2,o = LiflLQfliljQfl, generically one has -^2,0 = (0(1), 0(1), 0(1/A), 0(1))"^ at O2. 

By definition, the normahzed vectors are ipi^t = fi,t'ipi,t- Using a cychc permutation, 
we find that V^2fc,o = (0(A), 0(A), 1 + 0(A), 0(A)y and' that /2,o(A) = 0(A) at O2. Using 
the permutation argument again, we derive that /i+2,o(-^)//i,o(-^) = 0(A) at O2 for even i. 
Therefore, one has f2k,o{X) = O(A^) at O2. Now the required multiphcities for the vector 
4'2k,o at O2 follow. Furthermore, since 4'2k+i,o = -^2fc'^2fc,07 one can check that generically 



/2fc+i,o(A)//2fc,o(A) = 0(1) and ^p2k^ 



1,0 



(0(1), 0(A), 0(1), 0(1))^ at the point O2. This 



establishes also the multiplicities for the vector ip2k+i,Q at O2. 

Having proved the proposition for t = 0, one can prove it for t > by using the formula 
ipi^t+i = Pi,t'^i,t- Note that it suffices to study the cases t = and t = 1 only. Consider, 
for example, the multiplicity of the function V'j,!,! at the point O2. Since ipi^i^i = {ipi,o,2 — 
V'i,o,4)/(a^i(l + Vi + Zi+i)), one can check that the multiplicity of the right-hand side at O2 
is 1 — for i = 2k and it is equal to —k for i = 2k + 1, i.e., V^j,!,! and ipifl^i have the same 
multiplicities at O2. Other cases are treated in a similar way. □ 
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